Nonlinear evolution of a morphological instability in a strained epitaxial film 
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A strained epitaxial film deposited on a deformable substrate undergoes a morphological instability 
relaxing the elastic energy by surface diffusion. The nonlinear and nonlocal dynamical equations of 
such films with wetting interactions are derived and solved numerically in two and three dimensions. 
Above some critical thickness, the surface evolves towards an array of islands separated by a wetting 
layer. The island chemical potential decreases with its volume, so that the system experiences a 
non-interrupted coarsening described by power laws with a marked dimension dependence. 

PACS numbers: 68.55.-a, 81.15.Aa, 68.35. Ct 



The dynamics of semiconductor thin films is under ac- 
tive scrutiny due to its importance for both fundamen- 
tal science and technological applications [TJ [2] . Indeed, 
thin film elastic instabilities lead to the self-organization 
of nanostructures [3] potentially useful e.g. for quan- 
tum dots, wires and electronic devices with specific con- 
finement properties [4]. A notorious experimental ex- 
ample is Si/Ge films on a Si substrate which exhibit 
a variety of structures such as pre-pyramids, pyramids, 
domes and huts [5j [6]. Such epitaxial films experience 
an elastic stress due to the misfit with the substrate 
which is relaxed by a morphological instability similar 
to the Asaro-Tiller-Grinfeld thermodynamical instability 
in solid-liquid interfaces [7] . This instability was first ob- 
served in experiments in helium at low temperature [8] 
and more generally in various solid interfaces [9j [TUl E] • 

Although the evolution of epitaxial films involves many 
complex phenomena regarding surface energy, intermix- 
ing and kinetic processes, we focus here on the main ef- 
fects ruling the dynamics of the morphological instability 
in strained films. The dynamics is ruled here by surface 
diffusion driven by the interplay between isotropic sur- 
face energy and elastic energy [12] [13] . When the film is 
infinitely thick or when the substrate is infinitely rigid, 
different theoretical [TJJ [15] and numerical p~6j [T71 [18j [T9] 
approaches revealed finite-time singularities enforced by 
elastic stress concentration which account for experi- 
ments in thick films [HJ [9] where dislocations can finally 
develop. However, these models can not describe exper- 
iments of thin films in the Stranski-Krastanov type of 
growth [5j [6] where the surface organizes smoothly into 
islands separated by a wetting layer and evolving with a 
coarsening dynamics under annealing [6] . A crucial issue 
for these systems is the wetting of the substrate by the 
film [20j [21] which is a good candidate for regularizing 
the dynamics of the instability. Indeed, crack singulari- 
ties were circumvented near the instability threshold by 
considering slope dependent wetting effects [22]. How- 
ever, the interplay between elastic relaxation, surface en- 
ergy and wetting interactions is still under active study 
[23] [24] and the description of the long term dynamics of 
the morphological instability in a thin strained film is an 



open issue. In this Letter, we present a model based on 
continuum elasticity which we solve numerically reveal- 
ing the existence of a non-interrupted island coarsening. 

We consider specifically a three dimensional (3D) dis- 
location free film deposited on a substrate with slightly 
different lattice parameters and with a priori different 
isotropic elastic properties. During annealing, the film 
shape h(x,y, t) changes by surface diffusion (no external 
flux nor evaporation). The boundary at z = h(x,y 1 t) is 
free while the film-substrate interface at z = is coherent. 
In the reference state, the film is flat and the elastic en- 
ergy density equals S° = Ef (df — a s ) 2 / a 2 s (l — Vf), where 
a a , E a and v a are the lattice parameter, Young modulus 
and Poisson ratio of the solid a, with a = f for the film, 
and s for the substrate. The dynamical equation of the 
film shape is then, see [T2| [T3]. 
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with D, a constant related to surface diffusion, and 
Vs, the surface gradient. Both elastic JF el and surface 
jpmrf _ j dr7(/i)^l + \Vh\ 2 free energies contribute to 
the surface chemical potential \i = J(^ rel +^ rsurf ) / Sh which 
reads 

(i = £[h] + l{h)K(h) + 7'W/Vl + |V/i| 2 , (2) 

with £[h], the elastic energy density computed at z = 
h(x,y,t), 7, the isotropic surface energy and ft, the free 
surface mean curvature. To account for wetting p~4j [22] 
[23] [24] , the surface energy 7 is supposed to be a function 
of the film height extrapolating from the bulk value 7/ 
when h—>oo, to some upper value when h—>0. Here, the 
wetting interactions are described by the characteristic 
length 6 and the strength c w >0, and the surface energy 
is written as j(h) =7/ [1 + c w f (h/S)] with some function 
/ going to zero at infinity. In the following we will use 
/(£) =exp(— £) when a specific form is needed. Finally, 
we set the length unit to be lo = S° /jf, the characteristic 
length of the instability with the corresponding time unit 
t = l%/D V . 

To compute elastic energies, we use the isotropic 
continuum framework with stresses <j" proportional to 
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Finally, in two dimensions (2D), Eq. Q reduces to 
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with summation over repeated indices, n,p,q = x,y,z, 
5 pq , the Kronecker symbol, and e^ q = ^(d q Up + d p u q ) — 
r] a S pq (Spi-\-Sp2) where u is the displacement with respect 
to the reference state commensurate with the substrate 
so that r]f = a? /a s — 1 and r] s = 0. Since the system is 
at mechanical equilibrium, it satisfies d q (jp q = with the 
following boundary conditions: u s — > when z — > — oo 
and is continuous at z = 0, whereas is continuous at 
at the free surface z = h(x, y, t) with 
n. To solve for elasticity, we use 
the thin film approximation [23] assuming the thickness 
h to be an order e smaller than the characteristic length 
Iq. Hence, in the film, we consider the rescaled variables 
Z = z/e and get u considering an expansion up to e 3 , u = 
Xln=o e n u( n )(x, Z). In the plane substrate however, 
elasticity is solved as usually using Fourier transforms 
with respect to r = {x,y}, J-[h] = (27r) -2 J dre lhr 'ft(r) . 
Eventually, we calculate the elastic energy £ — \e^ q (7^ q 
up to e 2 , the first nonlinear term, and obtain the central 
equation of this Letter describing the film dynamics, 
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with i, j, = x : y. In Q, we use the notation Oijij = 

1? @xxyy @yyxx @xyyx @yxxy and 

^ijfcz = otherwise, and define the functionals Hij[h] = 
^{(kikj/tyFlh]} with k = |k|. The different elas- 
tic constants are uj\ = 2Ef(l — v 2 )/E s (l — v$\ 002 

n 



2v s ){l + v s )/E s {l-v f ) and 
UJ2 = 2E 2 (l - v 2 ) 2 /E 2 (l - v f ) 2 {l + u f ), which match 
2(1 + v eq ) in the case of equal film and substrate elastic 
properties, v s = Vf = v eq and E s =Ef = E eq . In the latter 
case, Eq. Q coincides with the linear analysis of [23]. In 
fact, up to order h 2 , Eq. Q can be derived thanks to the 
elastic free energy 
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where x-indices denote ^-derivatives and H is the 
Hilbert transform acting in Fourier space as TL[h x ] = 
{\k\J-[h]}. In the case of equal film and substrate 
elastic properties and without wetting, we retrieve the 
result of [19] describing a 2D semi-infinite film. 

We now investigate the dynamics predicted by Eqs. Q 
and (|6|. In the linear regime, considering small pertur- 
bations of amplitude exp[<j(k)t + zk-r] around a flat film 
of height ho, we find <r(k) = —ah 2 + uj\k z — bk 4 with 
a = 1 + c w f(ho/S) and b = c w f"(ho/5)/5 2 . Hence, when 
/"(£) is decreasing and positive, there exists some crit- 
ical height h c below which cr(k) < everywhere so that 
the film is linearly stable thanks to the wetting inter- 
actions. However, for ho > h c , <r(k) > in a given k- 
interval, and the film is linearly unstable. For small wet- 
ting length S and exponential wetting potential, one gets 
h c ~ — C02S hi(uoi5 2 / Ac w ). When c w =0, the film is always 
linearly unstable and the nonlinear numerical computa- 
tions exhibit generic finite-time blow-up solutions [22] . 








FIG. 1: Space-time evolution of a 2D film according to (|6| 
with ho = 0.1. Surface diffusion induces a non- interrupted 
coarsening until only one island is left surrounded by a wetting 
layer with height h wl . 

To analyze the nonlinear evolution in the presence 
of wetting effects, we performed numerical simulations 
using a pseudo-spectral method in a periodic box of 
length L. To be specific, we selected parameters depict- 
ing a Sio.sGeo.2 film on a Si substrate with Vf = 0.278, 
cji=2.44, 002 = 2.52 and 0^=2.34, leading to / = 200nm 
and with the value of the diffusion parameter D given 
in [13], to = 8 hours at 750°C, see [6 . In fact, thanks 




to space and time rescaling, only Vf and uo^/^z are rel- 
evant for characterizing Eqs. Q and The wetting 
potential is described in an indicative way by c w = 0.05 
and 5 = 0.005. The initial condition is a flat film per- 
turbed by a small noise with a mean initial height fto- 
As shown in Figs. [I] and [2j a film with fto > h c is first 
destabilized by the morphological instability which gen- 
erates surface undulations according to the linear growth. 
The linear stage is then quickly replaced by a nonlinear 
one characterized by the emergence of well-defined is- 
lands which grow without moving and with a decreasing 
aspect ratio. The islands are surrounded by a wetting 
layer smaller than h c which allows surface transport and 
the subsequent islands ripening. Hence, no singularity 
appears here in the nonlinear dynamics of the wetting 
film. Studying parity in h of the different terms in ([5]), 
one concludes following [19] that the last nonlinear non- 
local term drives the surface towards smooth peaks and 
deepening and sharpening valleys which would lead to 
singularities [19]. Here, the wetting effects included in 
7 (ft) enforce a higher energetic cost for small ft and thus 
stabilize the thin film. In fact, both nonlocal nonlinearity 
and wetting are needed to regularize the dynamics of the 
instability which we now characterize by its final state 
and time dependence. Note that a steady island dynam- 
ics was also found in [24] which however did not account 
consistently for the difference in the film and substrate 
elasticity and lead to different long time behavior. 

Within the present model, we observe that the system 
evolves continuously towards an equilibrium state char- 
acterized when fto > ft c by a single stable island in equi- 
librium with a wetting layer of height h wl in both 2D and 
3D, see e.g. Fig. [I] whereas when fto <ft c , the final stage 
is a flat film of height h wl = ho. The equilibrium prop- 
erties {h w \ island volume V measured above h w \ etc.) 
depend only on the homogeneous chemical potential ([2| 
and on the sign of fto — ft c for large enough L. Computing 
H and V as parametric functions of the film volume V^, 
we find that when ft<ft c , ^ = 7' (ft) increases with un- 
til Vj = L 2 h c , whereas when ft > ft c , \i depends only on V 
and is monotonously decreasing in both 2D and 3D, see 
Fig. [3] Hence, in a regime of well-separated islands, big- 



ger ones should always grow by surface diffusion at the 
expense of smaller ones. At equilibrium, we also com- 
pute the maximum height ft ma;E as function of the initial 
height fto, see Fig. [4] The system undergoes a discon- 
tinuous bifurcation as the difference h max — ho displays 
a jump at the transition height ft c which agrees within 
a few percents with the linear estimate ft c ~ 0.036 cor- 
responding to 7nm. This first-order like transition also 
shown in the fJ>(V) plot of Fig. [3] is at stake in similar 
instabilities [251 . 
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FIG. 3: Equilibrium surface chemical potential ([2]) as function 
of the flat film volume V = V^ when ft<ft c , and of the island 
volume V when h>h c in 2D (a), 3D (b). 




0.03 0.04 0.05 0.03 01)35 0.04 



FIG. 4: Maximal height of an equilibrium island as function 
of the initial height in 2D (a), 3D (b). 

Finally, to describe the dynamics of the island growth, 
we compute the surface roughness w(t) = [(h 2 ) — (ft) 2 ] 1 / 2 
and number of islands N(t). Both 2D and 3D simulations 
reveal a non-interrupted coarsening with power-law be- 
havior w(t)~tP and N(t) ~ l/t c , see Figs. [| and [6) For 
2D systems, we find /3 = 0.26 and £ = 0.59 over nearly 
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three decades. Similarly, over the last time-decade of the 
3D simulations, we find (3 = 1.3 and ( = 2.0 which are 
noticeably departing from the 2D values, illustrating the 
difference between diffusion process over a one or two di- 
mensional surface. If one sought a self-similar solution of 
Eqs. @ or ^ as h^t? H(r /t 1 ^) for large t, one would 
get (3 = 1/3, thence requiring a more elaborate theory ac- 
counting for the wetting layer and nonlocal nonlinearity. 
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FIG. 5: Roughness as function of time with L = 6700 in 2D 
(a) and L=104 in 3D (b). 
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FIG. 6: Evolution of the number of islands as in Fig. [5] 

In summary, we derived nonlinear and nonlocal equa- 
tions describing the stress driven morphological instabil- 
ity of a thin film on a deformable substrate with a priori 
different elastic constants and which account for wetting 
interactions. When both nonlocal nonlinearity and wet- 
ting are present, numerical simulations reveal a steady 
evolution towards an equilibrium state contrarily to the 
crack solutions predicted for the bulk morphological in- 
stability. When the film initial height is higher than some 
critical value given by the wetting interactions, the final 
stage consists of a single island with a chemical poten- 
tial monotonously decreasing with its volume. Consis- 
tently, the system undergoes a non-interrupted coarsen- 
ing in both two and three dimensions characterized by a 
power-law decrease of the island number with time which 
strongly depends on the system dimensionality. Fur- 
ther experiments on the number of islands of annealing 
films in the prepyramid regime of the Stranski-Krastanov 
mode [6] would be of great interest. New effects such as 
anisotropy and faceting will be included in future work. 
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